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- This peper ia concerned. with the following probl sn: 

Given @ p -dimenslonal normi random variable with means 
sere, yariancer one, and. correlation matrix KF 3; truncate 

thie random varieble in oll coordinstes, say at 1, tors t, 
respectively, end find expressions for E(%,"%,") after 
trumsation. An explicsi? seluticon of this spies is obteined 
for m=zi,2, nwd and ==i,n21, thet is for the 
exnectatiions. varisnocs and acvaricnaes of the distribution 
after trunsation. and an extension of ths method for greater 
values of m,n {a indicated. 


le Introduction. 
In varicua fielde of appifed atetiatias, aush aa naycho-= 
- logical meseuremer:ts and personnel seiectian, one frequently 


as sm ah wm A ~~ 2 med wom DT oe 
deals sith popeist sienzs which may Ss comsicsrsa os originaiiy 


al 
1s. Peer ae 
= PoSeentsd to the Tnstiture of Kathematical SUACLATLSE, 


Pbicago, Devexber 23, 1956. 


= nesearsh dens unter tne apensorsiiip os bes OFTice oi 
Havel Research, — 
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mpltiveriats normal, but cocdisied by truncation in each coordie 
nate seperately. Por example, @ p timsnsicnai normal paridon 
variable xX, 9 Ko pease ay may represent p quantitative 
traits of an individuals very often an admission teat requires 
thet each of these traits se ebove a certain presassiined 
Vaiue, so that only thoze individuais pass the twst for whom 
X, >t, seoo, X, 2, - It has been show [1], that this 
method of ‘saaebtine ten sous undesirable properties; it 4s 
howaver frequently the only practical method, and hence it | 
may be of come interest to study the properties cf digtributions 
obteined by such truncation, 

in the present paper explicit expressions are obtained 
fer the moments B(X,) » E(X,") , B(X,X,) » and it is indtegted 
how the method can be extended to the general case of = 
(Z,° 2,7) . Te possitility of truncation in sem but not 
sli soorfinetcs 1s ineluésd since e.g. the cass of X, not 


truneated, X, truncated at T corresponds to +t, = -co » 


ae ~~) s ).ae ee at 


i € Bepliess expressions ers slse cbtsined for the 
marginal probability density function of X, and for the 
joint marginal pedeZ. of (%, » Xp) » after truncation in 

- Exreunier are siven for the use of acm cf 
“tha results for deteraining t, , ty oeoc-. a0 that certain 
preeassigued wnanges in tne population ars achieved. | 


2, A ikmown lemma en detemuinests. . 
let R bea pxp matrix with the elemnmts ¥, 3 3 lst 
B,, be the eofantor of Py > RY the metrs= of the Bay ~ 


= 3 bead 
Hs 5 . the cofastar of Ry q ia R* 9 end oe Fghig? 
{(p~2) dfimeneicnal minor in n® obtained by deleting the 


i eth aad j «th rows and « -th and v <th solumms. Then 


we hare 


(21) fa™] = fap 
(2.2) May F a jnjP* 
(263) WE gage GRITS Cog tay oF yy Fay) IBIP? 


Sis proof cf this lemme may be found in stendaré treatises 
74 éeterninants, GeEo iS] Pedle 


3. Eqextions. for the moxents BX" x) - 
We al @ miitiener=s:1 ; -dimensional ranion varicble 


Zz, ry Sa aessea — e wate. tas ‘Sosrciaticn imivran RS iF, Pu ° 
asf @ iy Zeseoe P » SDE (without loss of generality) ae meee 
0 and variances 1. Itep.d.f.is 5 


(Sol) FR epee ney) cote oN 


The distribution is assumed to be non~singuler and nence the 
quadratic form > > + = ZX, 2 poadtive definite. 

a > / 
fruneasing ba My » Xe pereg x at ty 2 Be peoos ty ie 
wespestively es have Tac ths new p.d.f. of Eye Fg soeee | 
aftae truncation . | 7 
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(See) Bik; oXaresehy) ms 
0 eluawhers. 
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With these astatieas (3.2) becoses 


(3.3) B(Ry eX geve00K,) SS 


GN ey otgeceost aR) FU Sze 202%, ) fee % 5 tye Ml gk scccsds 


iv) Si BONNETS.» 
—— as Wf a a 
Sa oe a ee -) Bus apnae WounouD 
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Integrating the ight side bir parts, we obtain from (3.4) 
after soma si:aplificstion ; ‘ 


Be 


Pp 4 e , 
< a . n = Pr | seM@le Ty 
(3.5) Ry B(x, "X5") +> Ry HN es RIE + | giles 


2 a Ry ty". | 
al = oO 
ve. (ty pecret, 5h) in)” 3 & fet / eee 
{ 


” arya oh ai on 


eee MARE APE RY li PRE eS RRNA ARI a TL A AS om 4 en te Yh ek © i a (ne nes mR en Set en en 


= Se 
SPAY EL EITC TE Parra es PT RT AE ED EES I ED 


fo evaluate the integral on the ripht side we apply the 


transformation 
= F,.% 

(3.63 v= a Ve : ue Biv 
is 


and obtain, using the iemna of yectian 3, 


(‘os ses -3)) A 43%"; 
- ; - 
wf RR eg tyt es BR av. a, 
Spun | 
fr 


R.. “Sap gat 
£5.8) 2, , = ze Pie Hid = La sm2s3 ,ccear 


\j. -eoe ae  e . 
Tha matrix ¢ lat a ia positive definite sinee (Ba) 


sa Was acauisd positive davinite. 


Replacing the subscripts 1,2 by ¥,8 respectively we 
obtain 


(369) Ry Bix ™X,") + > art - IRIE ty (Xp X47) = 


(lS, @ Feyty) - 7 q 

oa : 

y “2 pw ; | 
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soe Bp 2Bl pZsevegP r i 

= 2 : 

(3-92) aye = aie Vie hy) (dare, iy imle2yscenD . 4 
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Tt oan be verified that the teverse of the watriz (2177) as | 

=e " 

(3092) Fi ike am Jai The} t. > ‘for Ly Iny2s seer : 1 
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that 4s the matrix of partial ocrrelation coefficients 


Ps gow e ae . a 


he Special cases: mel.,anwzGirme2,nnu02mnene21. 


rettsne m=ei,ne d in (3.9) ws obtein 


a (t.)! ¥ 5 
(hel) > Baw %! m Se Eo J nis. 
bel 
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We can rewrite the cquations (4.2) move conetsels, as 


ie ete -vcatL asatt ta |itnc-cher hchleshammmmiaiie lieu] Feel A emeah, Liicsanadaanimaseleadaslennsmlinindscieaste Uiccadeacca LL ete rr ee ree iateciaeaeae eee daceetl a tease cee 


To solve this ayste= of equations we uze (2.1) and (2.2) and 


obtain 
P 
Cio) By) = 2 Paw BiL) 5 WEly2roreoD 6 


ext, to obtein equations for the second moments, ws set 
me2,n20,end se#il,ns1 in (3.9). sing the notation 


win) 0 oft) Ver} caer 
: 3 ena ey : St 2 - vee 
G,(ty reo. ot 5R) (2 w) oles . f 
é 
: age t 
23 T Al) vv, 
iui, ya "Tbe ad ad os 
ste d @ -rey %y + Fyuty) ° us *e 
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we obtain, respootivelys 


(505) > RegBiz Ze) = ii BYRLIW) » Wel 2r0cesP 5 and 
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Combining (4.6) end (ye?) » we. have 


> 
(heB} $7 Rew(X A) = GRE Cae + Blwe)] , sel e2s0000p 
322 . mel eeeewe se e 


‘fhis ine aystex of 3° equations in the pipes UNnKNOwS — 


BRA) = BE) - Sines p~ 2, aig) , 


eefsis? a 


(none « mibersven of an gee independent equations, ‘The 
Sioas for hich iacsewa«yp form sucha systems to 


 qhow this, we arrencse these equations and their unimovwns ia tis 
manner indicated by the following table: | : 


it will be sufficient’ 


a 


eA 23 
w 


oO. 


Biz.2.) 
SO ee ee 

hy eg te 
hs 


OBO OF SSSE HOOHOESESSSSOSHSOSTE SOO TODOS OHSS SSEESCOSTEPHSHSHOSOSBOOHODOOHSE 
pasts : : see Pete ees : - ae 
OOo dere ccceeeseeseoesseoeuseeseTeSeSss i oseeseesecessreceessenoece 
sare : aa aes : aa 
BORSSOROSHSEE POETS OO RAFODEE AHS UH HORSE ZOSUTS = 


o se the he 


_ SFestuwETT ~~ 
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.In (4.9) ali the colums exoert the last sontein the 


eaeffiolents of the unknown indicated in the solum heeding, 
While the lest cvolum (headed £E) contains the richt side terms 
| or tnose equations uf system (455) for which ji SSsP®SPpe 


The determinant of the coefficient matrix in (4.9) 18 


Rr. 
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RS 
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7 ee oe Be 
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° 
e 
eee a 
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and is #0 since each factor is a principal miner of a 
poeitive definite matrix. Thus these syquations yield a solte 


tion for (4.8). Using (2.1) and (2.2) it 18 easily verified that 


= “~ a i= ae ot 


; BUR hy i= “fe yg itt 39) + fay 9 J9unn 


ia evaluating h(ws) we can use the previous resulta, 


in particular (4.4), on the firsé tere 


sissies Gia wld ae a hewn {el pp aetna al eeRchay rsad e iS 


Sa cs GND Retina ala anh aR al AR ae tal RU Rt abet cata « pa a Aa SL a + 
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Ce he SM Ne MME eT re en 1 i Rea gs eNO EE wae Si AE eh egies Savy te, Ee Oe ee Ee 


Ps » 12 = 
ree =o 1 P 
= ee wk Vg¥ 
(ese) » J 
eee ] \ ae v, e ' ya T av, ® 
433 ae ww #e 
Lore, 


for it represents, exeept for appropriate constant, the 
marginal expectation of the coordinate Y, of a (pel) 
déimensionel rendom variable Vy0 Voseces Vy we,?. Vv we19eee? Vy P 


truncated at ists : f2"F2uty wotes bye"? wed )ww 


> £ 2 
vw \f-r ae Vee (ek }w 
tr t, : 
eh (erliw ow , epee » Yespectively; the second fern 
oe bots 8 7 Tp 


_ oi prs 


G.(t 3R) (2m 2 t-Fyaty 
yoooest,, 9 


foo 4,3= j 
Tee se aw 
eee "sa's @ pn avs 
ter ¢ ole 


ott. = 


ye ter  ¢ i 
co i ae > 2\. 
& ip! *y 9eo0nt, 38) may, pS Al 1. / 


Bxpresaicus (ih) eee appear to be useful for - 
getting up and solving various practical probleme of the kind | 
iliustrated in Section 6. The numerical evaluation of these 
SRpTessions Ssquirss ths Gompiitation of integrals of the typs 
\Seaaje The Values of such iztegrais for pe2 may be found 
in Pearson's Table V1li. ~ 1X in (61, For pm3 and pel a largo 
tuumber of the required integrals may st{11 t< founda in theses 
tables, For ml ell such integrals involved have to be 
ceicv.iated, a cask which may require the use of high-speed 
 Gomputing equipment. 
Po obtain valuos of higher moments, ons sast go bask to 
£5.89 end by si:ilar manipulations as above, obtain the required 
nusiber of independent equations ¢o solve for the unknoms, | 


5e Expressions for the nargine) distributions of X, and 
(X,eX5) efter trunestion in XL, Xopeccs X,- 

zr WM) is the ped.f. of X, after truneation in 
“ae 2 acces Rs then a (3.3) 


F on | 
iS) seeert,, 3B) ee FE, vibe gXgiGkige» Gk, o es 
"| te adh , : 


Using the transformation 


gee Se =n emerasts dined iiemmanssenytalhedenmons ami santa ygichantndaeaippiclanasenconiiateatees easmpamereeer steer taterers 
Fic oct einen hese ont eee EE ER TORTIE 


a [ 
£.8 
sa Vy = cheek, sn 3 
; aa a segisccesP e¢2 | 
: | 
one obtains | 
RVs) =zxC, Jm2p3yeceeP 
4 de ncw ®) ail 8 Jer s3e0005P 
(503) 
‘ 44 - P31" 13 
BAT 5) = s Ly F223 ,ecesP 
if Lee, )(Lar% 4) 
vi B(x, v3? = 0 | Jm2a35-00sP e 
; a By (23, po313y Xy» Voocoos¥,, 22 again distributed 
according to the mlticaormal law, and hence according to 


(5.2) and (5.3), expression (5.3) beccsos: 


Ox. ) st e>.. x 4 


t_or .X. 
' G, (ty oocoet, sR) . \ fiers, yi-r?, . 


where Ty 4@ defined in (3.92). 


a of ERI By . If Yo (XoXo) denctes the p.d.f. of X, Xo after 
teuncation in Xj, Xozsees xy » then by (3.3) 


, : ie CO ; 
(t.5} GA.t%. 5X5} 4 ag 2 gk) [... J £( Ky soe) GRgeotX, @ 
= 3 ty 


§ 
Eis 


| 
| 
| 
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teing the transformation  ~ : 
Hey 
Xp = Xe 
Vg = rN TO 4% +Ou% + rs oo » I8SelpooesP » 
ag" & : 
ners - : 
"wy ial . | 
4, a Pox » ie oe _ 8g 123 slheeeesP 
ia 
eni a. is’ the cofector of r,, in 4 g » one easily verifies 
thas i ; 
By?) oo 2 ae 
(5.7) E(V.X,) a O ' Le3pleccoeP 
BAY sXo) =0 
Riv.v.) Sa. r re 
ech Salt ig i or way waar wear © & 55743 % zis +4 e Pod s 
¥ M11 43 O53 45 
i : ds Je3oocseyDe 


Henee, since (XysRos¥groeee¥y } again has e multivariate 
normal atrieuticn, we obéein from (546) and 5.7) 


77 Zar dnd tes Satis 


ir oa 5 x “G eS — 9 


UaltyeTyoe) = on iA, t, + + A,.x, + A aide ® Le3 pheesesP 
JO,,0, a ees oe 


Se (9,5) = (H(VAV,)) - 


b. Some Applications. 

The following probiem is of practical interest: ‘for a 
bivariate normal rantom variable (Xj XQ) with expectations o Py 
variances 1 and known correlation coefficient r , it is 

required to find t, end t, so that, after truncation at 
t, end t, , the expectations of XZ, and Z, assume the ; 
préewassigned values m and Ma 
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To find such %, » to » wo have according to (kok) 


{6.2 = Bil) + v &(2) 
Fh an — ae —_ BSD 
Vea) iti = aye TF Bas ° 


 Veing expression (eas for h(i) and simplifying, one obtains: 


: 
z oS eis. 4 


$$$ es eS Ee a EY: SS 
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aod ? 
(6.3) fas sr OE of two curves in the t, » tp plane, and 
“wila serva to illustrate ths procedure 


 $puncation ™ = 1.5 9 Mp = 2.0. The rightehand sides in 


SE SES TOT A FATTER et OO RCI 


; e & ¥ Vie? F ti. « 
6k} a: = EK = , : 


fhese equations show that the inequalities ro, < Mm, < § 2 mo, 
are neseesary for the existence of a solution. 
Ra jonny rumerical Values for t, , ty, » one may consider 


determine their intorssetion, Ths following mumerical exampie 


oe 


Given x = .60 and the required expectations after 


(6.3) and (6.4) become, respsctively, 3.67 and 0.469. By 
trial, using tebles, one finds the following three points on 
each of the curves: | 


ty ltt) «=, ty Rah ty) 
[ +35 | 1.532 | 3-67 | [.so] 2.430 | .469 


60 | 1.594 | 3.67 | 1.530 | .169 
wos | 1.655 | 3.67 | 


stot | es | 


Plotting these values one finds for the point of intersection 
t, = -554 » to = 16537 , and substituting these values into (6.1) 
and (6.2) one obtains m.* = 1,501 , m," = 2,001 which ta a 
good approximation to the vequirsd values. ; 

Wext we consider the following problem: We wish to truncats 
i and X, at t, and t, respectively so that the expectation 
of X,. after truncation has & pre-assigued value m, and the 


retained part of the population 


ee Ee nS ee ee ee 
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2 . 2 2. 
owe (X>° 2p x¥+ ¥) 
»2(Lex*} dx dy 


Git sto3 7) « sie Te ee 
a ewes 


ia ss large as possible. This is equivalent to maximising 
B(tysto3r) under the condition (6.1). Using Lagrange multipliers, 


we sonsiéger 


Ea A a TARTS A) AM NNT GENT 


B to+rt 
H(t, ,t,) = G(t, stor) + Alm,G(t..tosr) - glt,) sa os - 


=rt 
=P Hts) ea ; 


we wisk to solve the equéetions 


It is easily verified that these equations become, respectively, 


(6-5) A= + 
ty - m 


ert. 
(6 63 ale Pe eee 
@ (domed) aM - © 
heo® 


os 


) toert ert 
(6.7) m, Gltystosr) = @(t,) O(a) - v H(t.) da) =O. 
4 At, st, elt, ti P(t, Nit 


«19% 
e i h 


From (6.5) and (6.6) we obtain: 


It is well kmown [3], [4], that for finite Uso, 
> yen > U ., From (6.8) we see, therefors, that our problem 
} 
has He solution with t, and t, both finite. The following 
four possibilities remain: 


aj r>G, ts -0 , t,< Go 
b) r>0, t, < CO, tp * =0 
6) re0,t, =-m, ty < o 
@) r<e0,% < © » ty = -0 ° 


2 it.) 
In the casus a) and «s) (4.7) yielda: rab ie} 


Plt) 


“while 4n esses b) and @) (6.7) becomes —=—-se = My e Sinse 


P(t, ) 
aes 48 a monotonic function of U , inoweasing from 0 to 
v 
ry e we reach the follewing conclusions: 
For r> 06, (cases a) end b) ) , max G(t, »tosr) under 
condition B(X,7) = & is obtained by truncating in x, alone 


Paap saree a aS IE CPs RRR PARSE ROR ESSERE MSIE BORE RAE SEES AE 1 ORR Ne EE SARL LISE AEA RIES SEN USNC IE BLE ER SERENADE TTT pelea e? 


«= 20 « 


- (t, ) 
at ty , wnerse ty is obteined fron ‘ = Ml) Only 
; t 
a hase | 


positive vaiuss of m can be achieved. 
For r<0O, (cases c) and dad) ), we mst truncata fn 
P(t.) 
X alons and use the solution obtained from ee ss a ° 
for m, <0 ; and trumcats in X, aeione, at t, obtained 
) aou) 
im weaereiedite: ME for m > 0. Tables of SAU) may be 
ye d(T) 
i ; 
found 6.g6 in i6]. 
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